Granular Brownian Motors: role of gas anisotropy and inelasticity 
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We investigate the motion of a 2D wedge-shaped object (a granular Brownian motor), which is 
restricted to move along the x-axis and cannot rotate, as gas particles collide with it. We show 
that its steady-state drift, resulting from inelastic gas-motor collisions, is dramatically affected by 
anisotropy in the velocity distribution of the gas. We identify the dimensionless parameter providing 
the dependence of this drift on shape, masses, inelasticity, and anisotropy: the anisotropy leads to 
dramatically enhanced drift of the motor, which should easily be visible in experimental realizations. 
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Introduction — We investigate the motion of a 2D 
wedge-shaped object, which we shall refer to as the mo- 
tor (Fig. [T]) . It cannot rotate and is restricted to move 
along the x-axis, as gas particles collide with it. When 
the motor experiences elastic collisions, there is a finite 
transient drift as the motor approaches thermal equilib- 
rium with the gas pQ. A finite steady-state motion is 
achieved when the gas- motor collisions are inelastic [5HB] . 
The latter systems have consequently been called granu- 
lar Brownian motors. 

They are prototypes of systems where small particles 
collide with heavy objects that break reflection symme- 
try. Such models have been used to explore the recti- 
fication of thermal fluctuations [21 [5], the adiabatic 
piston [21 [TO], and have lead to a novel treatment of 
non-equilibrium steady-states In an experimental 

realization [5], it was demonstrated that they even obey 
non-equilibrium fluctuation theorems. 

So far, however, all pertinent theoretical studies are 
based on thermostatted gasses such that impacting par- 
ticles are sampled from a Maxwellian velocity distribu- 
tion. When thermostatting via stochastic forcing, this is 
a reasonable assumption [3]. On the other hand, exper- 
imental realizations of granular gasses typically exhibit 
sustained heterogeneities in density and granular tem- 
perature [12H16] . Moreover, when shaking in the plane 
of observation, they exhibit noticeable anisotropy of the 
granular temperature [17j . Consequentially we denote 
them as anisotropic gasses. 

Here, we revisit the approach by which [H [7J [S] derived 
the theory for the isotropic case. Then we address the 
motion of the motor driven by an anisotropic gas. 

Gas Velocity Distribution Function (VDF) — Follow- 
ing |17j . we model an anisotropic VDF using a squeezed 
Gaussian, 
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where m is the particle mass, k is Boltzmann's constant, 
T = (m(v x + Vy)), /(2/c) is the gas temperature aver- 
aged over both degrees of freedom, and T x and T y are 
the granular temperatures in the x and y direction, re- 
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FIG. 1. A particle (black circle) colliding with the Brownian 
motor (triangular wedge with wedge angle 2#o). The angles of 
the edges, i G {0, 1,2}, are measured counter-clockwise from 
the positive a;-axis to the outside of the motor, yielding 9o, 
6\ — tv — 8q, and Bi = 3n/2, respectively. 



spectively. Anisotropy is quantified via the squeezing pa- 
rameter, a 2 := T y /T x . a > 1 as we only address vertical 
shaking. 

Introducing dimensionless velocities, v :— v/ yj2kT '/to, 
and requiring cj>(v x , v y ) dv x dv y = 4> a (v x ,v y )dv x dv y , re- 
duces the VDF to, 
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which depends on a only, and not on m, k, T x and T y . 

Gas-Particle Interaction — A collision event is illus- 
trated in Fig. [I] The motor has dimensionless velocity 
V = Ve x , and a mass M. Collision rules depend on 
which side of the motor, i g {0, 1, 2}, is being impacted 
and on the coefficient of restitution, r. 

Assuming no change in the tangential component of 
the gas particles velocity, 
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where tj = (cos 8i, sin is the tangential vector to the 
surface being impacted. In contrast, due to restitution 
the reflection law for the normal direction becomes, 
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where hi = (sin 6i, — cos Qi) is the normal vector. Single 
collisions obey conservation of momentum, 

v' x + MV = v x +MV , (2c) 

where Ai := M/m is the mass ratio. Altogether Eqs. ^ 
determine the change in the motor velocity, 

Ui := V - V = 7 l (v x - V - v v cot Qi) , (3a) 

where 

7, = y(r,MA) ■= (1 + r) ,ff °\ „ ■ (3b) 

M + sin 0j 

Time Evolution of the Motor VDF — For indepen- 
dent collisions, the probability density, Pt(V), of finding 
a motor with velocity V at time t, follows the master 



equation, 

d t P t (V)= [ W(V -u;u)P t {V -u)du 

- / W(V;-u)P t (V)du, (4) 

where W(V; u) du is the conditional probability of a mo- 
tor experiencing a collision resulting in a velocity change 
V — > V + u. It can be expressed as an integral involving 
four specifications: selecting only those outcomes which 
are (i) commensurate with single collisions (Eqs. (|3|), 
and (ii) collide with the outside of the motor's surface; 
(iii) weighting single particle collisions by the impact fre- 
quency, where the collision frequency for a stationary mo- 
tor is used to non-dimensionalize time; (iv) sampling over 
all possible impact speeds and the motor's sides, where 
Wi(6o) is the probability of picking the side i [2]: 
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Consequentially the steady-state solutions of Eq. Q are 
selected by a, j(r, M.,0), and the wedge angle 28q. 

Solutions to the Master Equation using Moment Hier- 
archies — Given that, Vn £ N + and m < n the deriva- 
tives d™(u n W(V;u)) vanish for u — > ±oo, the Kramers- 
Moyal |18j expansion can be applied to the moments, 
M fe (t) :='(V k ) = J R V k P{V 7 t)dV. Together with the 
jump moments, a n (V) :— L u n W(V; u) du, we arrive at 
an evolution equation for the moments: 
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In order to accommodate a more general velocity dis- 
tribution, we compute the jump moments by expanding 
them as a power series, 



i=0 

such that Eq. ^ reduces to an infinite linear system, 

oo 
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reminiscent of a matrix equation with matrix elements, 
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Time-resolved motor velocity-PDF — In general one 
still can not solve the infinite matrix equation Hence, 
we truncate Eq. ^ at order N, which leads to, 
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The expansion coefficients d n ^ in Eq. (|7|) are computed 



using the Taylor expansion coefficients d r , 

,(<) 



where an' (V) is the z-th derivative of the n-th jump mo- 
ment. In order to compute these derivatives, the delta- 
distribution in Eq. ([5| is integrated out, resulting in non- 
trivial integrals. As long as V = 0, these can be evaluated 
using Mathematica. The higher order derivatives of these 
integrals are related to each other allowing them to be 
computed recursively. This provides an analytical, albeit 
tedious, expression for Eq. (10 1. 

Asymptotic analysis reveals that d n .i ~ — i~ 1 ^ 2 for 
large z, resulting in a combined truncation error in 
Eq. (fl0|) of the order of 10" 10 for N = 20. In this work, 



we hence solve Eq. ( 10 1 for N = 20 and a wedge angle 
$o = 7r /4 unless stated otherwise. The initial condition 
will always be an ensemble where all the motors are at 
rest: M(0) = (1,0,0, •••)• 

Fig. [2] illustrates typical time dependencies of the 
motor drift, (V), and motor temperature, T '■= 

M (j^V 2 ) — {V) 2 ^j . (i) For elastic collisions and an 
isotropic gas, the ensemble undergoes a finite transient 
drift while it heats up to the temperature of the gas [T|. 
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We observe that, for small T, 




FIG. 2. (Color online) Ensemble drift, (V) (main panel), 
and temperature, T (inset), against t for motors with mass 
ratio A4 = 10, and 8q — n/4. Blue lines (i): r = 1 and 
a = 1, elastic collisions with an isotropic gas. Green lines 
(it): r — 0.3 and a = 1, strongly inelastic collisions with an 
isotropic gas. The motor relaxes to the values predicted by [2] 
(black horizontal lines). Red lines (Hi): r = 1.0 and a — 1.02, 
elastic collisions with a slightly anisotropic gas. 



Subsequently, the drift ceases, (ii) When introducing in- 
elastic gas-motor collisions, the steady-state acquires a 
finite drift velocity and a temperature significantly lower 
than the gas [2J. (Hi) Here we note that a small amount 
of squeezing, a — 1.02, causes a drift similar to the drift 
in a system with strongly inelastic collisions. Note that 
this squeezing hardly affects the temperature. 

In the subsequent sections, we examine the parame- 
ter dependence of the steady-state drift, (V), and motor 
temperature, T, respectively. 

Motor Drift — The inset in Fig. [3] (a) shows that for 
a fixed coefficient of restitution (r = 0.3), the drift ve- 
locity initially scales as 1/A4. For large M. and a / 1 
it approaches a constant value depending only on a and 
6>o- The l/OVf scaling is in agreement with the theory for 
the isotropic gas [2j. We conclude that the drift for light 
motors is affected primarily by the inelastic nature of the 
gas-motor interactions. Here the theory for the isotropic 
gas is a good approximation. In contrast, massive motors 
are more strongly influenced by the anisotropy of the gas, 
no matter how slight this may be. 

In order to fully characterize this crossover, we consider 
the limit of a massive motor: A4 — > oo. In this limit the 
7i term in Eq. (3a) simplifies, 
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Due to this factorization of sinfl and T, massive mo- 
tors undergoing dissipative collisions (r < 1) behave like 
motors undergoing elastic collisions (r = 1) yet with a 
slightly higher mass. This is in agreement with results 
for the granular Boltzmann equation pjJJ [20]. Conse- 
quentially, the limit of a massive motor corresponds to 
the limit T — > + and is independent of restitution, r. 



dn A ^° r 
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Hence, for isotropic gas VDFs (where a = 1), the matrix 
defined by Eq. Q becomes upper-triangular in leading 
order of T. This corresponds to the decoupling of the 
time-evolution equations for the moments, as observed in 
[2j. In contrast, for a > 1, the time evolution equations 
for the moments become coupled again: 
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This shall be the starting point of a perturbation the- 
ory around (T, a) = (0 + ,l). We assume that, in the 
limit r — > + the steady state is still largely indepen- 
dent of truncation size for small (a — 1). Hence, we find 
that the null space of the upper left 2x2 sub-matrix of 
Eq. ( 13 1 accurately determines the steady state drift due 



to anisotropy, 



di,o 
d%A 



(sin0o-l)(a-l). 



(14) 



Note that Eq. ( 14 ) does not depend on M. . This is quite 



astounding since it implies that the drift velocity of the 
massive motor is of the order of the gas-particle veloc- 
ity (dimensionless (^) ani is of the order 1), even though 
the transferred momentum from the gas remains constant 
with increasing Ai. 

The crossover occurs when the drift for the isotropic 
case (^)iso - (l-r-)X- 1 v /^/2(sin6'o-l)/4 is of the 
same order as the drift due to anisotropy. Consequently 
the dimensionless number, 
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(V) s 



(V)i 
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~r : 



(15) 



characterizes the dominant driving of the motor. For 
(3 <C 1, the dynamics is driven by inelastic collisions (r < 
1), and for (3 ^> 1 the dynamics is driven by anisotropy 
(a > 1). Plotting (V)/(V) an i as a function of /3 provides 
an excellent data collapse, Fig. [3] (a). 

Motor Temperature — Fig. [3] (b) shows that the tem- 
perature is independent of M. for M. > 10 and it is 
affected by inelastic collisions more severely than by 
anisotropy. We now follow the perturbation theory of 
the previous section to determine the correction to T in 
first order of (a — 1). 

Since the motor temperature contains a coefficient of 
1/r, we must expand A to second order in T. According 
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icr 2 io 1 io° 10 1 io 2 io 3 io 4 10 1 io 2 io 3 io 4 

(3 = 4M(a-l)/(l-r) l/T=M/(l+r) 

FIG. 3. (Color online) Data for all combinations of 8 = tt/4, r S {0.3,0.5,0.8}, a G {1.02, 1.007, 1.002, 1.0007, 1.0002}; and 
9o = 7T/10, r = 0.5, a G {1.02, 1.007, 1.002, 1.0007, 1.0002}. (a) Master plot for the motor drift where the inset illustrates of the 
effect of varying a. The curves for inelastic collisions with an isotropic gas: r — 0.3, a = 1 (straight blue line, (i)); and elastic 
collisions with an anisotropic gas: r = 1, a = 1.02 (straight red line, (ii)) have been included for reference. All other curves 
show the drift for decreasing a G {1.02, 1.007, 1.002, 1.0007, 1.0002} from top to bottom, (b) (top) The motor temperature, T, 
for 0o = 7r/4, r G {0.3,0.5,0.8}, a G {1.02, 1.007, 1.002, 1.0007, 1.0002}. (bottom) The difference between motor temperature 
and the asymptotic theory. For comparison, 6 Q = tt/10, r = 0.5, a G {1.02, 1.007, 1.002, 1.0007, 1.0002} is also shown (★). 



to Eqs. (12 1 A then takes the form, 
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This results in a further increase of the coupling between 
the different moments. In order to reliably compute 
(V 2 ) an i, the null-space of at least the upper left 4x4 sub- 



matrix of Eq. ( 16 1 must be used, yielding the asymptotic 



1 



expression for the temperature, 



Tani^l+ ^p(l-sin(9 ) 2 + sin 2 6> (a - 1). 

(17) 

The lower panel of Fig. [3] (b), shows the converge onto 
this asymptotic value. 

Conclusion — We have investigated the motion of a 
granular Brownian motor that is driven by inelastic col- 
lisions (particle- motor coefficient of restitution r) with an 
anisotropic velocity distribution (with anisotropy a—1), 
modelled using a squeezed Gaussian, Eq. ([!]). 

Examining the scaling of the drift with relative motor 
mass, M, we identified a crossover from the motor drift 
arising due to inelastic gas-motor collisions, to a setting 
where it arises predominantly from the anisotropy of the 
gas. Examining the steady-state drift of the motor in 
the limit of large M. , we have identified a dimensionless 



parameter /3, Eq. (15) (independent of wedge angle). For 
ft <C 1 inelastic collisions drive the drift of the motor, and 
anisotropy is negligible; for /3 S> 1 anisotropy dominates 
the drift and restitution in motor-gas collisions becomes 
negligible. In the latter regime we have identified a re- 
markably strong enhancement of the drift: it is of the 
order of gas particle velocity, even in the limit of infinite 
motor-particle mass ratios. Is this remarkable regime ac- 
cessible experimentally? 

Many experiments, involving agitated granular mat- 
ter, are kept in a steady state via shaking from the 
walls. Such systems always exhibit an anisotropic veloc- 
ity distribution j!7j . Laboratory experiments can have an 
anisotropy of the order of a « 2 [3T] , and the most conser- 
vative estimate for simulations yields a ~ 1.12 ([T7j Fig 
4, inset). Given maximally inelastic collisions (r close 
to 0) this amounts to (3 ~ 0.5M. For M > 10 typical 
experimental realizations therefore probe, at best, the 
crossover regime rather than a regime where the drift 
solely arises from the inelastic collisions. If one wishes 
to probe the latter regime, isotropy of the gas particles 
must be enhanced by at least two orders of magnitude 
for the experimental setups we are aware of. 

The dramatic enhancement of the drift thus lies in an 
easily accessible regime, and it certainly calls for further 
experimental and numerical exploration. 

We are grateful to P. Colberg, S. Herminghaus, R. 
Kapral, W. Losert, D. van der Meer, L. Rondoni, and 
M. Schroter for enlightening discussions. 
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